H 2 dissociation over Au-nanowires and the fractional conductance quantum 
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The dissociation of H2 molecules over stretched Au nanowires and its effect on the conductance 
are analyzed using a combination of Density Functional (DFT) total energy calculations and non- 
equilibrium Keldysh-Green function methods. Our DFT simulations reproduce the characteristic 
formation of Au monoatomic chains with a conductance close to Go — 2e 2 /h. These stretched Au 
nanowires are shown to be better catalysts for H2 dissociation than Au surfaces. This is confirmed 
by the nanowire conductance evidence: while insensitive to molecular hydrogen, atomic hydrogen 
induces the appearance of fractional conductances (G ~ 0.5Go) as observed experimentally. 
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Gold surfaces are chemically inert and regarded as poor 
catalysts at variance with other metal surfaces. The low 
reactivity of molecular hydrogen on noble metal surfaces, 
like Au and Cu, seems to be well understood 0,0,0,3- 
Density Functional (DFT) calculations Q have shown 
that the dissociation of H2 on Au or Cu is an activated 
process: for large molecule-surface distance, d, the in- 
teraction energy is repulsive with a high barrier of 1.1 
eV (around d = 1.5 A), which the molecule has to over- 
come to move along the reaction path H2 — > H + H and 
reach the final atomic chemisorption state, with a to- 
tal adsorption energy (2.07 eV per atom that is less 
than the H2 binding energy ( 4.75 eV). Compared to 
surfaces, small particles are known to be better cata- 
lysts. Au is, in particular, considered as an exceptional 
catalyst when prepared as nanoparticles on a variety of 
support materials 5J. Understanding this strong catalytic 
activity is still the subject of an extensive research effort 
with different possible explanations, including the par- 
ticle shape or p^imeter^ support effects and the metal 
oxidation state 



Nanowires are a good example of systems whose size 
is so small that one can expect their reactivity with 
molecules to be increased considerably. The forma- 
tion of metallic nanocontacts has been analyzed in de- 
tail thanks to the gentle control of the distance at the 
atomic scale provided by both the Scanning Tunnel- 
ing Microscope and the Mecanically Controllable Break 
Junction LU|. While, in many metallic contacts, the 
formation of nanowires during the last stages of the 
stretching, just before the breaking point, is character- 
ized by an atomic dimer geometry (Al is a paradigmatic 
case [HI El ) •> in Au the final geometry seems to be a 
chain of several atoms between the two electrodes with 
a conductance close to the conductance quantum unit 
G = ^f- 0. Although the formation and stability of 
monoatomic Au chains has been addressed by several au- 
thors 0,0,0? certain relevant aspects -in particular, 



the changes in the structural and transport properties 
of the nanocontacts induced by the presence of impuri- 
ties HE BE El- are not yet fully understood. 

Recently , Csonka et al |2l[ have analyzed the inter- 
action of H2 with a breaking gold nanowire and have 
found new fractional peaks (in units of G Q ) in the con- 
ductance histogram. Moreover, conductance traces in 
a stretched nanowire demonstrate a reversible transition 
between fractional and integer conductances, in a time 
scale of milliseconds or seconds, suggesting successive ad- 
sorption and desorption of hydrogen on the chain. These 
experiments do not show these effects for Cu and Ag, 
where stable single-atom chains are not formed. This 
suggests the great importance of the Au chains in the 
variation of gold conductance in the presence of H2. 
Csonka et al 21] have proposed a possible explanation 
for the observed behaviour in terms of a dimerization ef- 
fect which has been theoretically predicted for idealized 
clean Au nanowires 0, [2^. However this dimeriza- 
tion has neither been observed in conductance histogram 
measurements nor in more complex theoretical simula- 
tions II EI El. 



In this Letter, we show that there is, in fact, a strong 
link between the enhanced reactivity of the stretched 
monoatomic gold chains and the appearance of the frac- 
tional conductance peaks. First, we have simulated the 
whole deformation process for an Au nanocontact upon 
stretching, finding the formation and the final breaking 
of a 4-atom Au chain. The realistic nanocontact config- 
urations, calculated in this way, are then used to investi- 
gate whether the new fractional peaks in the conductance 
are associated with adsorbed molecular or atomic hydro- 
gen. As our analysis suggests that only atomic hydro- 
gen can be responsible for these changes in the nanowire 
conductance, we have also investigated how the chemi- 
cal reaction H2 — > H + H is affected by the presence of 
a freely suspended Au-wire. This simplified model cap- 
tures the key ingredients in the real nanocontact struc- 
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FIG. 1: Total energy per unit cell and ball-and-stick structure 
models (see also Figure 2) for the Au nanowire as a function 
of the stretching displacement. The inset shows the total 
differential conductance and channel contribution along the 
stretching path. 

ture and provides a natural playground to explore the 
influence of the low dimensionality and the strain in the 
different steps of the dissociation process. These calcu- 
lations show that a stretched Au-nanowire is much more 
reactive than the Au surfaces, with a small activation 
barrier, around 0.1 eV, for the H2 dissociation and larger 
chemisorption energies. Our results for the nanocontact 
conductance, combined with the low value we have cal- 
culated for the i^-reaction activation barrier, strongly 
suggest that the molecule dissociates on a Au-nanowire 
and that the observed fractional conductance upon ad- 
sorption of molecular hydrogen is basically due to the 
atomic hydrogen produced in the reaction. 

Our calculations for stretched Au-nanowires have been 
performed using a fast local-orbital DFT-LDA code 
(Fireball2004 24]). This code offers a very favorable 
accuracy /efficiency balance if the basis of excited pseu- 
doatomic orbitals 25] is chosen carefully. The electrical 
conductance of the Au nanocontacts has been calculated 
using a Keldysh- Green function approach based on the 
first-principles tight-binding Hamiltonian obtained from 
the Fireball code, at each point of the deformation path 
(see references [ll|,|l2l for details). 

First, we have analyzed the formation of a Au- 
nanowire obtained by stretching a thick Au wire hav- 
ing four layers, with three atoms in each layer, sand- 
wiched between two (lll)-oriented metal electrodes as 
shown in figure ^ (configuration A represents the ini- 
tial relaxed configuration). We use a supercell approach, 
where periodic boundary conditions along different di- 
rections are introduced: parallel to the surface we have 
considered a 3x3-periodicity, while in the perpendicular 
direction we join artificially the last layers of both elec- 




FIG. 2: Ball-and-stick structure models and total conduc- 
tance for the adsorption of molecular (top) or atomic hydro- 
gen (bottom) on different nanowire geometries (corresponding 
to the points E,F and G in Figure Q}. 

trodes(see refs. 0,^2)- Figure shows the total energy 
of the system as a function of the stretching displace- 
ment 26]: notice the energy jumps associated with the 
irreversible deformations and structural rearrangements 
the wire has during this process. This figure (see also 
fig. 2) shows several snapshots for different geometries 
corresponding to the labels in the energy curve. Our 
DFT simulations reproduce the characteristic formation 
of Au monoatomic chains (with up to four atoms) found 
in the experiments. Similar results have been recently 
obtained using parametrized tight-binding molecular dy- 
namics |27| . 

The inset of figure [I] shows the conductance of the sys- 
tem along the stretching process, as well as the differ- 
ent channels contributing to it. Notice that these results 
compare well with the experimental evidence . In par- 
ticular, we reproduce (a) the long conductance plateau 
associated with the formation of the monoatomic chain 
(configurations C — > I), where the conductance is ba- 
sically controlled by a single channel associated mostly 
with the Au s-electrons; and (b) the conductance oscil- 
lations during the elongation process. The very good 
agreement between both our structural and conductance 
results for the evolution of the Au nanocontact and the 
experimental evidence provides strong support to the re- 
maining simulations presented in this paper. 

In a second step, we have analyzed the nanowire con- 
ductance upon the adsorption of molecular and atomic 
hydrogen. Starting with the different geometries corre- 
sponding to the points E, F and G in figure d we have 
analyzed, via our local-orbital DFT code, how molecular 
and atomic hydrogen are adsorbed on those geometries 
and, then, how the nanowire conductance is modified ac- 
cording to the new optimized structures (see Figure [2j. 
Molecular hydrogen is weakly adsorbed on all the dif- 
ferent nanowires, with energies of around 0.3 eV. Notice 
that, in all these cases, molecular hydrogen does not pen- 
etrate the nanowire too much (see figure EJ) and the con- 
ductance properties of the nanowire are not affected prac- 
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FIG. 3: (a) Density of states (DOS) for H chemisorbed on the 
Au(lll) surface (solid line) and on the Au nanowire (dashed 
line), (b) Comparison of the DOS for topmost atoms in the 
Au(lll) surface and for Au atoms in the four atom chain of 
the nanowire (see inset). 

tically by the adsorption of molecular hydrogen. In par- 
ticular, for the three cases shown in figureEl the nanowire 
conductance takes the values 1.05 G G , 1.05 G Q and 0.95 
G , respectively. Atomic hydrogen introduces more dra- 
matic changes: in particular, for the cases shown in fig- 
ure [2j the adsorption energies are 3.3 eV ,3.5 eV and 
3.9 eV for the chains with two, three or four atoms in 
the nanowire, respectively. These values are much larger 
than the ones found for H adsorbed on a surface (around 
2.1 eV) 0,11. Moreover, we find significant modifications 
in the conductance of these three cases, with total val- 
ues of 0.68 G , 0.65 G and 0.55 G G , respectively. The 
eigenchannel analysis for all these cases shows that the 
transport is dominated, as in the clean nanowire, by a 
channel mostly associated with the Au s-electrons but 
with a reduced transmitivity. This reduction is related 
to the significant displacement of the density of states 
(DOS) to lower energies (particularly evident for the d 
bands) for the Au atoms bonded to hydrogen that results 
in a reduction of the total DOS at Ep. The nanowires 
with three or four atoms have been reanalyzed assuming 
that two hydrogen atoms are simultaneously adsorbed 
on the chain: the case of three atoms (case F) presents 
a conductance of 0.45 G G , while for a four- atom nanowire 
(case G) the conductance is 0.2 G Q . 

The enhanced reactivity of Au-chains, with respect to 
Au-surfaces, is due to the change in the Au DOS. Fig. 3b 
shows this DOS for the atoms of the Au(lll)-surface, and 
for the Au- atoms of the four atom chain in the nanowire: 
Au-atoms with lower coordination form directional bonds 
and present a narrower DOS shifted towards the Fermi 
level. Then, for H chemisorbed on the chain, its DOS (see 
fig. 3a) presents a tightly bound state at -7.3 eV below 
the Fermi level, whereas the DOS for H on the Au(lll)- 
surface presents a broadened resonance with some con- 
tribution from antibonding states just below the Fermi 
level J]| . Notice that the reactivity of Au chains is further 




FIG. 4: Energetics of the H2 dissociation on a stretched Au 
nanowire along an arbitrary configurational path. 

increased by the wire stretching due to a further shift of 
the Au-bands towards the Fermi level. 

It is remarkable that, in the geometries shown in fig- 
ure [3 neither the molecule nor the atom penetrates the 
nanowire, breaking the bond between two Au atoms. Re- 
cent work by Barnett et al 17] showed that, for an essen- 
tially broken Au wire (with G ~ 0.02Go), a barrierless 
insertion of the H2 molecule into the contact is possi- 
ble. This would correspond, in our case, to a wire state 
beyond the configuration I in figure Q where the con- 
ductance is still very close to G Q . They then used this 
configuration as a starting point for a detailed study of 
the structure, orientation and stability of the molecule 
upon compression of the wire. Our approach is different, 
since we are interested in the wire-H2 interaction during 
the chain formation process. We have calculated the en- 
ergy barriers the molecule experiences when moving from 
the geometries shown in figure El to inserted sites. Our 
calculations yield values larger than 0.5 eV. 

Both the conductance and the total energy results dis- 
cussed so far support the dissociation of molecular hy- 
drogen in the Au monoatomic chains. As the energetics 
were discussed in terms of the local orbital code, using 
the LDA approximation for the exchange-correlation, we 
have re-analyzed the dissociation mechanism on freely 
suspended A u- wires, with H (or H2) adsorbed on them, 
using CASTEP j^l H^, with a gradient corrected 
approximation [^} (GGA) for the exchange-correlation 
functional. This simplified model for the nanowire is dic- 
tated by the computational resources needed for this full 
calculation, but it offers the possibility to discuss the 
relative contribution of the low dimensionality and the 
strain in the enhanced reactivity of the nanowire. 

We have considered a freely suspended Au wire with 
six independent atoms per chain and periodic boundary 
conditions chosen to produce a stretch deformation in or- 
der to simulate the breaking process. After relaxing the 
free wire for each strain condition, we have studied how 
molecular hydrogen interacts with it: in particular, we 
have explored the possibility of having the following re- 
action, H2 — > H + H on the suspended wire. Figure 0] 
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shows our main results for a wire with a 11% strain in 
an arbitrary configurational path: the case A in the fig- 
ure represents the geometry of H2 interacting with the 
freely suspended wire; this is a physisorbed state with 
an adsorption energy of around 0.02 eV. Case B corre- 
sponds to an intermediate state in which a H-atom is 
chemisorbed between two Au-atoms , and the other H 
is still bonded to the Au-atom on which the molecule 
was initially physisorbed: this case has an adsorption en- 
ergy of 0.3 eV; we have found, however, a barrier of 0.1 
eV between states A and B . Case C corresponds to the 
final reaction state, in which the two H are adsorbed be- 
tween two Au-atoms; the chemisorption energy of this 
final state is 1.4 eV, but again we find an energy bar- 
rier of 0.4 eV between states B and C. The reaction path 
drawn in figure 0] shows that molecular hydrogen sees 
a total barrier of 0.1 eV for the reaction H2 — > H + H on 
a freely suspended Au-wire. In a Au-surface, the energy 
barrier for that reaction is around 1.0 eV, and the adsorp- 
tion energy, with respect to H2 , is negative, around 0.5 
eV. These numbers show the extreme importance that a 
stretched Au-nanowire has in the reactivity of molecular 
hydrogen: its reaction energy barrier and its chemisor- 
tion energy are lowered at least by 1 eV by the nanowire. 
Notice that the reduced dimensionality of the nanowire 
(compared to the surface) is not enough to induced this 
high reactivity, as shown by similar calculations for a 
non-strained chain, where we find an initial bound state 
for the molecule with energy -0.10 eV and quite large 
barriers (~ 0.5 eV) for the dissociation. 

In conclusion, stretched Au nanowires are much more 
reactive with molecular hydrogen than Au surfaces. Our 
DFT-GGA calculations show that, in a stretched Au 
nanowire, the activation barrier for H2 dissociation is 
very small, around 0.10 eV. Notice that this is a reli- 
able upper bound, as the use of the the GGA approx- 
imation (that only partially corrects the gross overesti- 
mation of barriers in LDA), the freely suspended wire 
geometry considered, and the neglect of possible quan- 
tum tunneling effects in H tend to overestimate the cal- 
culated barrier. This enhanced reactivity can also be ex- 
pected in the case of Pt, where the formation of chains of 
several atoms has been also observed l32|. Complemen- 
tary evidence is given by Csonka et al |2l[ data for the 
conductance of a stretched Au nanowire with adsorbed 
molecular hydrogen. Our calculations indicate that only 
atomic hydrogen can be responsible of the changes ob- 
served by those researchers in the nanowire conductance. 
It is the combination of these two results, our calculated 
activation barrier and the changes in the nanowire con- 
ductance, that strongly suggest that molecular hydrogen 
dissociates when adsorbed on stretched Au nanowires. 
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